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ABSTRACT 


The normal impingement of plane turbulent jets on smooth walls 
is studied experimentally and analytically. Three distinct flow regions 
exist in this phenomenon; the free jet region, the impingement region 
and the wall jeu reguonm, Lie free jer region exhibits properties aimost 
Tdeneicalelo these Oh anicee eG. Lis extent 1s show to be ew fixed 
portion of the impingement height. The impingement region is characterized 
by increaséd static pressure and continuous deflection of the velocity 
VEClor. wNelvelocity ana pressure fields an this region are studied 
by means of a Pitot cylinder. A measuring technique, not involving 
rotation and separate pressure measurement, is developed. Analytical 
expressions are developed for the variation of axial velocity and 
pressure as well as for the maximum velocity along the wall and the wall 
Shear stress, the flew in the wall get region approximates that. of the 
classical wall jet. The length scale varies linearly with distance. The 
local skin friction coefficient is shown to be a constant, whereas the 
velocity scale is’ shown to vary as a power of distance, the exponent 
being less than -0.5. Generally, the dimensionless height of impingement 
has little significance, while the nozzle Reynolds number is a relevant 


parameter. 
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CHAPTER I 
INTRODUCTION 


The impingement of a turbulent jet on a solid surface finds 
application in a number of engineering problems. Jets issuing from 
hydraulic outlet works, weirs, vertical take-off aircraft and various 
spraying devices are examples of such problems. 

The very general case of the impingement problem is a jet 
issuing from an outlet of any shape and impinging on a solid surface of 
any shape and texture. However, plane and circular jets impinging on 
plane surfaces may be recognized to be, by far, of more practical im- 
portance. Furthermore, since the plane jet impingement constitutes the 
simplest case of this prone it is obvious that its study should be 
of importance in understanding more complicated situations. 

This thesis is an attempt to study both experimentally and 
analytically, the case of a plane turbulent submerged jet impinging 
normally On’a smooth and flat boundary: ‘In this case, the flow issuing 
from a two-dimensional nozzle strikes a flat, smooth surface as shown 
in Fig. 1.1. After impingement, the jet changes direction and flows 
along the boundary. 

In general, we may recognize three distinct regions of flow 
Wien reference to its pattern |Lj, (2), [315 Im region 1; “the7ziiow 


characteristics are, for all practical purposes, identical to those of 
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the "free jet region". In region II, the flow undergoes deflection and 


at the end of the’ region becomes parallel to the wall (region III), 
assuming a pattern very similar to what is known as the wall jet. It 
is appropriate therefore to refer to region II as the "impingement 
region" and to region III sf the “wall jet region". It is logical to 
expect transitions between these regions. This aspect will be dis- 
cussed in later chapters. 

The main objective for region I was to determine its extent. 
In region II, where the static pressure is generally higher than the 
hydrostatic value and the velocity vector assumes a relatively large 


angle with the initial jet direction, the main objective was to study 
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the velocity and pressure fields. The wall shear stress was also 
studied. In the wall jet region the objectives were to investigate 
the velocity field, the jet thickness and the distribution of the wall 
shear stress. 

A number of investigations regarding impinging jets have 
been carried out already and some experimental results as well as ana- 
byetGaleprocedures ale avaliable. (U), (215) [oasis tlD lo miiece are 
mostly connected with circular jets and include extended studies of 
the free jet and wall jet regions as well as some measurements of the 
axial velocity in the impingement region [2]. For the plane case, 
Mathieu [5] reported measurements for impingement at small angles. 

His work is mainly concerned with the wall jet forming after impinge- 
ment. Also for the plane case, Schauer and Eustis [1] reported re- 
sults for different angles of impingement including the case of normal 
impingement. Their work is also concerned mainly with the wall jet 
region. The wall pressure and the maximum velocity,in the impingement 
region, along the wall were measured. 

The present investigation purports to provide detailed re- 
sults for the impingement region where very little is known and to en- 
large the scope of results already available in the wall jet and free 


qee regions. 
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CHAPTER IT 
EXPERIMENTS 


Ze Experimental Set-Up 

The experiments were conducted in the arrangement shown in 
Pigs. 2.-and 2.2 usine air as the flow medium. 

The air was supplied by a compressor connected to a rec- 
tangular plenum chamber which terminated in a nozzle smoothly tapered 
EOvansexit width of 0.083 in. (Fie. 222).)) Ihe aspect ratio was about 
65 so that a plane jet could be produced. 

In the chamber, a series of circular tubes was placed near 
the entrance and a sequence of five screens was placed further down- 
stream, as shown in Fig. 2.2. The tubes served the purpose of 
"straightening" the flow streamlines which tended to become curvi- 
linear because of the sudden expansion at the entrance. The screens 
were placed in order to produce complete mixing of the flow before it 
entered the nozzle section of the chamber. The intensity of the flow 
was adjusted with the aid of a pressure regulatory valve. The 
velocity at the nozzle exit was varied between 115 and 207 fps. In 
this range, compressibility effects are known to be negligible. 

To produce impingement, a 1/4 in. thick plexiglass plate, 

6 in. wide and 5.75 ft. long was placed normal to the axial plane of 
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the jet. One foot high sidewalls were used to avoid spilling. The 
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5 
sidewalls were also made of 1/4 in. thick plexiglass and were of the 
same length as the plate (Fig. 2.1). The impingement height, i.e. 
the distance between the plate and the exit of the nozzle, was adjust- 
ed by raising or lowering the chamber and was varied between 0.103 
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2.2 Experimental Techniques 
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Referring to Fig. 2.3, the Bernoulli theorem applied between 


section 0-0 and 1-1 gives: 
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where 910 sis the mass density (of ithestludid. Py> UL and Py? U 


are the static pressure and velocity at the sections 1-1 and 0-0 
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The validity of Eq. 2.2 was checked by comparing manometer heights with 
nozzle velocity heads at different velocities Us - No difference was 
recorded between corresponding readings. Since the manometer was 
vertical and contained water, the arene reading would be about 3 in. 
Assuming that reading errors were at most +0.05 in., the maximum head 
recording error would be about 1.6%. Since he ) YP, » the maxi- 
mum error in obtaining ke values would be less than 1%. 

Thus, the nozzle velocity, Us » was measured indirectly by 
recording the pressure Ap, . This procedure has the advantage that 
Ue may be recorded within seconds at any instant, so that corrections 
for small supply fluctuations can be applied. Manometer heights were 
converted to pressure and velocity assuming a temperature of 20°C and 


an atmospheric pressure of 14.7 psi. 


2.2.2 Velocity and Pressure Measurements in the Impingement Region. 


Pitot Cylinder. 
In the impingement region the direction of the velocity vector 


is not known before hand. In addition, the static pressure is higher 
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than the hydrostatic. Therefore, use of ordinary velocity measuring 
devices» was) ruled out. ©After a séarch of relevant literature’ [6], [7], 


[8], [9], it was decided that a Pitot-cylinder be used. 
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(b) TOP VIEW 


The Pitot cylinder is a cylindrical instrument of small 
diameter having three orifices (see Fig. 2.4) near its edge which is 
rounded. When the instrument is placed in any flow and the orifices 
are connected to a manometer, three readings will be obtained. The 
values of these readings will depend upon the velocity, V , the static 
pressure, Po >; and the direction of the velocity vector with respect 
to the axis of the cylinder. When the flow is two-dimensional the 
instrument may be placed perpendicular to the plane of the flow. Then, 


the direction of the velocity vector will be fully described by the 
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8 
angle , as shown in Fig. 2.4. The ordinary procedure is to rotate 
the instrument until readings 1 and 3 become equal. Then, reading 2 
represents the total head and the direction of V is indicated. If 
the static pressure at the point is known, the velocity may be calcu- 
lated. To avoid rotation of the instrument and separate pressure 
measurement, an alternative procedure was developed, as follows; if 
Py» Po> P3 are the pressures recorded for each of the orifices 1, 2, 3 


respectively, we may in general write: 


Dee ee ei ee ee ee) (223) 


where 9 is the mass density of the fluid and K, are three cali- 
bration factors. Ordinarily, these factors will depend upon the value 
of @ and the tube Reynolds number, R = Vd/v (wv is the kine- 
matic viscosity). However, from what is known about pressure distri- 
bution around circular cylinders [10], [11], dependence of K, on 

Ro should be negligible for not excessively large values of 9. 


Thus, it is assumed that: 
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From Eqs. 2.3 we obtain: 
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where h,, ho and he are corresponding piezometric heads. Using 


Eqs. 2.3 1C 1s €asy to show that 
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of gd may be found from Eq. 2.7 by measuring hy; hos h, and form- 
h, - h, 
pRgvede (ale Veba pepe Ga eh) eee Knowing ¢$ , the values of K,, K,, K, may 
hy = h, Ne Gers 


be determined from (2.4). The values of V and Py may be calcu- 


lated from the following equations: 
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The above equations may easily be derived from Eqs. 2.3. To evaluate 
the functions fi = 1,2,3) the following procedure was employed. 

A Pitot cylinder 1/8 in. in diameter with orifices located at 1/2 in. 
from the edge,was calibrated in the wind tunnel of the Hydraulics 
Laboratory for three different values of Re . The angle between 
successive orifices was 45°. The velocity and the static pressure 
were held constant in each run, while the instrument was rotated by 
incremental angles and the piezometric heads, h., were recorded 
each time. 

The values of K. were computed as indicated in Eqs. 2.5, 
and are plotted against 6 in Fig. 2.5. This figure indicates that 
for the present range of Ro values (1430 to 5500), any dependence 
of K. on a is indeed negligible. The present results can be 
extended to larger or smaller values of Rp» when using geometrically 
similar instruments, for values of © so that | >| <0 30c es Lhas ere— 
striction is imposed to ensure that none of the orifices will be 
located in the wake zone where the Reynolds number has significant 
effect on the pressure distribution. A lower limit for R could be 


fi 


taken as 300, since above this value the wake is fully turbulent. An 
upper limit for Re is obviously the eriticalwelue about Jax 10° : 


however, it is doubtful that values greater than 50 to 60 thousand 


will be encountered in practice. 
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The curves of Fig. 2.5 are almost symmetrical about the 
axis $ = 0°. Deviations from symmetry are believed to be due to 
imperfections in manufacturing. The calibration factor K was com- 
puted from the above curves and is plotted against $ in Fig. 2.6. 
Figure 2.7 shows the variation of the factors Ky ~ Ky and Ky - K, 
which appear-an Eqs. 2.8 and 2.8A. 

POr che present study the: chants of Pies .s2.5,. (2.6 candi, 
were used for better accuracy, since the tube which was used for 
calibration purposes was also used for subsequent measurements. Use 
Of Ehenpresent, technique may be done in conjunction with Pies. 2.3 to 
2.10, where symmetry is taken into account and generalized curves are 
drawn. 

Equations 2.8 were found to agree very closely. It is 
recommended that both are used and an average value of ho be adopted. 
In cases where one of the factors K, - Ky and Ky = Ky has a very 
small absolute value, only the equation involving the other factor 
sHowld be vused. | Regarding has. .2.95 it is recommended that the one 
involving the highest absolute value of K, be used for better 
accuracy. The equation involving the next higher value of K, may 
also be used as a check. 

The technique described above was based on the same line of 
thought as that adopted by Rajaratnam and Muralidhar [9] for cali- 
brating yaw and pitch probes, and was mainly developed for measure- 


ments in the impingement region. It was also used for measurements 


in the free jet region where the velocity field is well known. 
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2 
Measurements in this region could be compared against available data 
as a means of checking out this technique. No correction for turbu- 
lence was applied. 

50 as to be able to insert the cylinder in the flow, slots 
were opened on the sidewalls of the flume as shown in Fig. Desh l8 he 
Roughness effects of the slots were overcome by taping the slots. 

The measurements were made in the neighborhood of the centerline of 
the flume. 

To measure the pressure distribution on the wall, 1/16 in. 
diameter, static taps, located as shown in Fig. 2.12, were connected 
to a precision manometer. Each distribution was measured two or three 
times because of the violent fluctuations of the pressure heads. 

Measurerents of maximum horizontal velocities at the edge 
of the boundary layer were made with the aid of a total head tube. 

The static pressure was assumed to be equal to the wall pressure at 
the same vertical. The tube had a rectangular opening with a height 
of 0.67 mm and a width of 3.7 mm. This tube was also used for measur- 
ing velocities in the wall jet where the static pressure was assumed 


to be atmospheric. 


2.2.3 Wall Shear Measurements 

Preston's technique [12] was employed to measure wall shear 
stress with the aid of the calibration curves given by Patel [13] for 
smooth boundaries and zero pressure gradient flows. This technique 


is briefly discussed below. 
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For a Preston tube of external diameter d Lesting on “4 
smooth boundary, if Ap is the difference between the total pressure 
indicated by the tube and the static pressure on the boundary, it 


could be shown that; 


oe oe ea) (2.10) 


where 9, V are respectively the mass density and kinematic viscosity 
Ole Chey fluid, a is the boundary shear stress and F denotes a 

: F 2 2 2 2 
functional relation. If we set Apd’/4pv = Ap, and td /4ov" = Tose 


we may write: 
be F(T Sx) (2211) 


Using the equations given by Patel for different ranges of the vari- 
ables, Ap, and Mere EGuwc.d tom lOtted sims creme lo. 

Ordinarily, this method is applicable when tubes of circular 
cross-sections are used. However, subsequent investigations [14], [15] 
have shown that flattened tubes can also be used if the dimensionless 
groups, Ap,, tee are formed using the total height of the tube for 
the external diameter. Since a flattened tube was used to measure 
horizontal velocities, it was decided that the reading, when the tube 


rested on the wall, be used for computing shear stress. To control 


this procedure, wall shear was also measured by means of a circular 
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14 
tube at several stations. This tube had an external diameter of 1.28 
mm and a ratio of internal to external diameter of 0.65. The shear 
velocities obtained thus were plotted against each other as shown in 
Fig. 2.14. This figure was used to convert shear velocity values 
obtained by the flattened tube into values that would be obtained by 
the circular tube. The difference between corresponding values is 
generally small. 

In the impingement region the pressure gradient is appreci- 
able causing therefore some uncertainty regarding shear measurements. 
Shear stress was measured by the above technique and subsequently 
Patel's criteria regarding reliability of the results were applied, 


as discussed in Chapter III. 


Wee Experimental Results 

Measurements were conducted for different nozzle velocities 
and different impingement heights. The nozzle width remained constant 
throughout the experiments. 

Table I shows the combinations of impingement heights and 
nozzle velocities used,as well as the type cf measurements performed 
in each case. The symbols appearing in this table are explained in 
the definition sketch, Fig. 2.15,as well as in the "List of Symbols”. 
Referring to Fig. 2.15, in order to retain conventional notations for 
the co-ordinate systems and the quantities involved, two different 
systems are used. The system x, y refers to the free jet and the 


impingement regions excluding the wail boundary layer, whereas the 
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16 
system x2 Vy refers to the wall jet region, including the boundary 
layer in the impingement region. Thus, x, xy and y, Vy always 
denote longitudinal and transverse distances respectively, while the 
letters u and v always denote longitudinal and transverse velocity 
components respectively. The suffix m generally denotes maximum at 
a particular cross-section. 

The experimental results, which follow, do not include 


finalized plots. These will be presented in the next chapter together 


with analytical considerations. 
Be ook Free Jet and Impingement Regions 


(i) The velocity field 
Detailed velocity measurements were made only in run No. l. 
Transverse distributions of the velocity component u are shown in 
Piece lO. err us is the axial velocity and Le is the value of y 
at which u = Su, , these profiles may be normalized by plotting 
u ey ¢ / ae : 
—— sas sas shown in Fie. 2.172 This figure sugeests that the 
m u 
profiles remain similar as far from the nozzle as 8572 of the impinge— 
femt Heleht, H, In this range, the error distribution: 
=(04834n)¢ 
£(n) e (al) 
(here "E(Q) = 


m 
least for n < 2.0. Tollmien's theoretical curve is also a good 
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Ay, 
description of the data although it is not shown. Due to considerable 
scatter occurring at small values of f(n), it is difficult to 
determine the value of mn for which f(n) = 0. However, the common- 
ly accepted value of n = 2.3 appears reasonable. Therefore, 

Be ole 

The profiles are not similar for = values greater than 
about 0.85, although the deviation is not large. It is interesting 
to note in Fig. 2.16 that the u components assume negative values 
at large values of y . As the wall is approached further, these com- 
ponents are reduced in magnitude and very near the wall they tend to 
disappear (x = 5.63 in.). This may be explained if it is understood 
that at large values of y , the velocity vector assumes the direction 
of the (time average) streamlines of the wall jet. Therefore, u_ be- 
comes negative and it is nothing but the Vy component of the wall 
jet velocity. .As the wall is approached, the angle of the streamline 
with the horizontal tends to zero and so does the magnitude of Nege 

Fig. 2.18 shows distributions of the horizontal velocity 
component, v . As would be expected, v acquires a maximum value 
and then decreases slowly to become the uy component of velocity in 
the wall jet. 

As the flow deflects in the impingement region, a boundary 
layer forms on the wall. The velocity vector is almost parallel to 


the wall and acquires a maximum value at the edge of the boundary 


layer. This quantity is denoted by ua and “a typical variation with 
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x is shown in Fig. 2.19. This resembles the curves of Fig. 2.18, 
although ud does not necessarily occur at a constant distance above 
the wall. 

the distribution ofithe product wusy “is shown in Fie. 2.20. 
Physically, uv if multiplied by p represents the flux of x-momentum 
in the y-direction. Figure 2.20 suggests that uv is almost independ- 
ent Or (xe for "y’<S.072 in. “As’we shalloseevin the next chapter, 
this observation is of considerable importance in predicting the vari- 
ation of ae with x in the impingement region. 

(ii) The Pressure Field 
Figure 2.21 shows pressure profiles in the impingement 


region. The symbol h stands for manometer reading at a slope 1 


10 
vertical to 10 horizontal. Reliable profiles for a 4S LOGAN) could 
not be obtained with the present equipment because the pressure was 
so small that the relative error became unacceptable. The term 
pressure will be used to denote the excess static pressure above the 
atmospheric in the remaining part of this thesis. If Pa is the 
axial pressure and ae is thesvalue of ~y Tat whzehs sp = =P, einen 


the pressure profiles may be reduced to a single curve if ee is 


plotted against ee This iseshowns in Piles, =2.22.) 216 os Le 
p m 
n= te the functions 
~(0.834n)- 
Aa eae (2.12A) 


describes the data reasonably well. 
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19 
figure 2,23, (shows a typical wall pressure distribution.) sin 
the next chapter it will be shown how all the wall pressure profiles 


may be reduced to a single curve. 


(iii) Length, Velocity and Pressure Scales 
We have seen already how the velocity and pressure profiles 
may be made similar. For full evaluation of the velocity and pressure 


fields the variation of the scales uw. be Dae b must be known. 
u 
Figure 2,24 shows the variation of ue with oe for 
fe) fo) 


different dimensionless heights The decay of the maximum 


2 
Oo 


velocity in the free jet is given by the equation: 


u Cy 
ao aoe (2213) 
m x - x! 
2b 
fe) 
where ¢, is a comstant with an average value of 2.35 and x' is a 


il 
virtual origin correction, usually equal to a few nozzle widths. The 
present data follow the above law for a certain distance,after which 
they deviate more and more from this law. The data are described by 
Bdge 2 oat Cy = 2 Oa eeicw =e. 3D (2b). ThereLore,. £orsthis eset 


of experiments: 


ig 2 
cone ai ae eres. (2.14) 
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This equation is applicable at distances greater than the length of 
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20 
the potential core which for the present data is about 8 nozzle widths. 
A typical distribution of the axial pressure, pes is 
shown in Fig. 2.25. The pressure is seen to increase rapidly for 
x 2 4.5 in., reaching a maximum value at the wall. 
Figures 2.26 and 2.27 show the variation of the length 
scales Li and ae with distance x . The scale dD is seen to 


Valve Linea iy wi Line cer Ore 6 xe ie following the equation: 
Da = 0.1(x + 0.9) @2e15) 


Equation 2.15 is almost identical with the free jet equation which 
states that be = 0209/7 (= + x"), where Keeler anvat bite tere re am 
correction, not necessarily eoecserencinte to the virtual origin for 
velocities. 
(iv) The Wall Shear Stress 
Figure 2.28 shows typical variations of shear velocity, 


along the wall. By definition u, = Vv T,/9 : 


with distance x * 


u,5 v 
The shear velocity, and therefore the shear stress, rises from zero 
at the axis to a maximum value, decreasing steadily beyond this point. 
Since these data were obtained by use of Patel's curves [13] for zero 
pressure gradient flows, their reliability is subject to certain 


limitations, because the pressure gradient is significant in this 


region. These limitations will be discussed in detail in the next 


chapter. 
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Zi Diee Wall Jet Region 
This region was investigated with respect to the velocity 
field, the wall shear stress and the growth of the jet thickness. 


Since the horizontal velocity component, is almost equal to the 


uy> 
total magnitude of the velocity vector and the vertical component, Vy ; 
is very small in this region, the velocity field was studied with 


respect to the u, component only. 
(i) Velocity Distribution 


Typical velocity distributions in the Or direction are 


shown, in Bie. 2.29, 


: , eal 
a Ich by is the value of V4 at which U) = > Ua and 
ae < o , the velocity profiles may be reduced to a single similarity 
ul u y 
curve if - is plotted against — 2 1bas Ss eshowmin ete. (2530; 
mL al 


Although, it is difficult to determine the value of ny at which 


£(n,) =. wa value n. = 1.8 appears to be reasonable. Thus 


nb 
Be = 1.8 b ~ the walue ot Ny at which £(n,) = 1 may be taken as 


ny = OF LO ens. 7) Oe). 6 by . The present results agree closely 
with those reported by Schauer and Eustis [1] . The distribution for 
the classical wall jet [16] is also shown and is seen to agree well 
with the data for 4 < 1.0. The discrepancy between the experimental 
data and the classical wall jet curve may be explained if it is recog— 
nized that the velocity uy is likely to assume negative values for 


ny > 1.8, which are nothing but entrainment velocities, wee for 


the initial jet. Such negative components are not present in the case 
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Ze 
of the classical wall jet, so that uy tends to zero at a somewhat 


larger value of "Nh According to the above, the velocity profile 


could be sketched as shown in Fig. 2.31. 
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Then, the "true" velocity scale would be Sees lvl and the "true" 


length scale would be somewhat larger than b Say) pel.) eine 
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VarLoacion. of at [vt with By should agree better with the 


classical wall jet curve. 

Measurements of Se components were not carried out in 
this work, although their existence was noticed. 

A transitional zone exists between the impingement and wall 
jet regions in which the flow undergoes a process towards establish- 
ment of similarity. It is difficult to define the distance after 


which the profiles become similar. This is because transition towards 
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22 
similarity takes place gradually. Careful inspection of the data 
suggests that the free=mixing region is established first (ite. at 


earlier xy stations),followed by establishment of similarity in the 


boundary layer. 
(ii) Velocity and Length Scales 


Figure 2.32 shows the decay of the dimensionless velocity 


u x 
ml , y : F i : 
scale SUL with dimensionless distance, ment from the stagnation 
oO me) 


7 


point. A dependence on the dimensionless height a sly elo nig Welt keys 
fo) 


Data reported by Schauer and Eustis [1], are also shown. These data 


were obtained at a value of nozzle Reynolds number of about 


R_ = 43,000. Since the present results were obtained for Ro less 
u 
than 10,000, a dependence of “ upon the Reynolds number, Ro? is 
fo) 
indicated. 


The variation of the length scale, by» with xy is shown 


in Fie. 2.33. The results indicate a linear variation. This is in 


agreement with what is known already about the wall jet [17, ch. 10]. 


The value of C. = “a is approximately 0.04 and appears to be in- 
1 
dependent of the parameter oe [hts conclusion is supported by, 
fe) 
H 
the results of [1], where C. Onn Om Lor oS 20 and 40. The 
O 


difference between the present value of C. siaicl jelavehin (epg. LIKI) alts 
possibly due to the large difference in the Reynolds number. An 


average value of C in classical wall jets may be taken as 0.068 
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fig rch. 10). 
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(iii) Wall Shear Stress 
This will be best studied in conjunction with analytical 


considerations in the next chapter. 
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CHAPTER III 
ANALYSIS 


This chapter is an attempt to produce a means of predicting 
different properties of the plane impinging jet with the aid of 
dimensional analysis and theoretical equations. Emphasis is placed on 
the study of the impingement region. Since a qualitative understand- 
ing of the problem precedes any analytical attempt, we propose to pre- 


sent briefly the various physical concepts involved. 


3.1 Physical Considerations 

The flow from the nozzle enters a region in which the main 
properties of the jet vary almost identically with those of the corres- 
ponding free jet. Thus, us varies as ae ; De Varies as) x ; 
Vv << u , and the velocity profiles are similar for x greater thau 
the length of the potential core. Therefore the effect of the wall 
in this region must be very small. This effect is transmitted by means 
of increased static pressures, hence the excess static pressures must 
be very small in comparison with dynamic pressures. 

Downstream of this region the presence of the wall signifi- 
cantly alters the above characteristics. The flow is deflected and 


becomes parallel to the wall, so that ultimately it forms a wall jet. 


Between the free jet and the wall jet we have the impingement region. 
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The impingement region is characterized by a sharp increase in the 
Static oe a fast decay of the vertical velocity component, u, 
and by the fact that the horizontal velocity component, v, is gen- 
erally of the same order of magnitude as u. Quantitative definition 
of the limits of this region is arbitrary because it depends upon the 
particular criteria used in defining those limits. Two particular 
boundaries are of distinct interest in this case. One is the be- 
ginning of the region, Xo» which is the boundary between the free 
jet and impingement regions at the axis. The other is the end of the 
region along the wall, which could be defined as that distance X19 
at which the pressure increase becomes very small. 

If the subscript f is used to denote properties of the 
corresponding free jet, then Xo could arbitrarily be defined as that 


Valle Of sa mabiwniche) Ul 10.90 Gu On bie =.102..b or 
m m u u 


f£ nf 


aes O05 eo The values of Xo defined thus may be expected to be 
close, they cannot, however, be expected to coincide. 

The impingement region is characterized by reduced turbulence 
near the axis. This may be deduced from the results of previous in- 
vestigations [2], 113], where 1b appears. that the total head) on the 
axis does not vary excessively between x = Xo ands xo= Ue, Lhe 
flow is thus almost potential in the neighborhood of the stagnation 
point [3] and a flow pattern known as "stagnation flow" [10, p. 87], 
is likely to occur. Under these conditions, the boundary layer could 


be laminar for some distance from the stagnation point. Further 
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ee 
downstream along the wall, the flow becomes established as a wall jet. 
A transition between the end of the impingement region and the wall 
jet is likely to exist. The main difference between the wall jet 
after impingement and the classical plane wall jet, lies in the manner 
in which they are produced. The latter is produced by the efflux of 
momentum through a small opening. This situation closely approximates 
the idealization of a point source. However, the wall jet after 
impingement is produced by a pressure distribution on the axis and a 
favourable pressure gradient persisting for a certain distance down- 


stream. 
3.2 Dimensional Analysis 


SRP aRe! Free Jet and Impingement Regions 
For a plane, turbulent free jet, it has been proved both 


analytically and experimentally that: 


UE? be = yg oe) (301) 


where M = 0-+2b aie = efflux of momentum. Equation 3.1 suggests 
O O 
that viscous effects are negligible and that the concept of a point 


SouLce: LS valid. 


Assuming that the above are valid in the case of an 


impinging jet also, we may write: 
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(322) 


(33) 


Applying the principles of dimensional analysis to Eqs. 3-2 and 3.3 


we obtain: 
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The function h@) is schematically sketched in Fig. 3.1. If the 


beginning of the impingement region is defined to occur at a distance 


=| x 
° 


x at which uy = 0.98 a then 1S) TENS Solhenzaiom Ore 


ie? 
6 
h() = 0.98. jTherefore, the distance Xo is jaetixed fraction or 
the impingement height, H. The experimental results are plotted in 
Fig. 3.2 which supports the arguments used in the dimensional analysis 
leading to Ba, 3.4. 
Eq. 3.6 may be applied to the axis (y = 0) to find the axial 
pressure, p. 
m H = x 
Gi 


Rea 
ou 7/2 2b | 3 


i ee) (3.9) 


The stagnation pressure Ps? will be; 
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Ps H 


a os 2. g, (1) = const. (3510) 
fo) 


Dividing (3.9) with (10) we obtain. 
a es 
m7 85 Gan 


The experimental observations on pressure variation are plotted in 
Pigs 3.9 and the structure of Eq. 3.11 is well supported. Althouch 
the scatter is appreciable it is believed to be due to the number of 
computational steps involved in obtaining values of pressure, as well 
as to the difficulty of obtaining De 

The wall pressure distribution, Pi.» may be obtained from 


Eq. 3.0 tor x =H. Then’: 


We tt pe 3 (1,4) = eG) (et) 


In Eq. 3.12, y° is the distance alone theswall trom the 
stagnation point. It may therefore be substituted by xy in the co- 
ordinate system Xo Vz ° Hence: 
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Dividing (3.13) with. (3.10) we obtain 
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Equation 3.14 is verified in Fig. 3.4. The data suggest: that the 
x 
pressure vanishes when x) = 0O.4H . However, even at — = 0.30, 


the pressure is only 5% of the stagnation pressure, Ps> and the 
gradient is very small. Hence, it is suggested that the end of the 


impingement region be taken as x Re Os H ae 


i 

We have seen already that the pressure profiles are similar 
in the impingement region. To check whether similarity extends to the 
wall we proceed as follows. The value of Me at, the wall, Dawe is 
EOUnd | irom the data Of Pag.) 3.4 to. bes 


The similarity cusve for pressures is) given by Eq. 112A: 


2. 
_ _ _-(0.834n) ae 
DUD ies SG)F se sf Wie te E 
Pp 
cal ail 
Ne iene wall, i = hs = 0.1340 = 7.45 (=)? . Thus, if the wall 
pw 


profile is Similar to the protiles in the jet, tne function: 


xy 2 
Pa (0.834 * 7.45(C—)) 
ee ee re H 
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should describe the data. 


After performing computations, this becomes: 
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(25) 


This is plotted in Fig. 3.4 and appears to be in good agreement with 
experimental results, hence supporting the “statement that similarity 
extends to the wall. 

MEq.73+5, “b  isvaneunspecified leneth scale. Glnrthe 


present case two length scales are involved, one for velocities, b , 
u 


and one for pressures, 2 = *Lhus 
re x Be x 
H 7 829 G@ and Ge = 8G a 


The experimental results are plotted in Figures 3.5 and 3.6, in the 
form suggested by Eqs. 3.16. Although there are no results for 
different heights and nozzle velocities to confirm (3.16), it is 
believed tharwvenese equations are valid because they are based on the 
same assumptions as Hqs. 3.4, 3.11 and 3.14 which were veritied, 

For the wall shear eee To: and the horizontal velocity 


component, u at the edge of the boundary layer in the impingement 
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region, we could also write: 
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The experimental data are plotted, as suggested in (nt7) on Figures 
3.7 and 3.8. Although there is some scatter, it is believed that Eqs. 
3.17 are valid with the possible exception of a Seealil region near the 
stagnation point where the boundary layer is likely to be laminar. 

In the impingement region the wall pressure gradient is 
appreciable, therefore, the reliability of shear stress data is sub- 
ject to limitations. 

Patel [13] suggests that in favourable pressure gradients, 
use of the Preston technique will produce errors not larger than 62 
if the following conditions are satisfied; 


GOO Ne 0.007 were |< 90-007 


u,d 
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Gia) wae ts 200 
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ae y 
(iii) Ax 0 
where A = — “2 » d = external tube diameter. These criteria 
pu, 


are rough guides rather than absolute restrictions. 

Condition (i) was derived empirically with the aid of fence 
readings. Condition (ii) is imposed to ensure that the tube is located 
in the region of wall similarity. Condition (iii) is imposed to en- 
sure that the flow is not subjected to relaminarization. In the 


present case, the only physically plausible sequence of different flow 
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62 
regimes is firstly laminar and later turbulent. If the boundary 
layer is laminar this must happen at the early stages of its develop- 
ment. Any sequence therefore, that includes laminar reversion in the 
direction of flow is physically unacceptable. 


Using Eqs.3.17B and 3.15 we find, after some manipulations, 


thats: 
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For condition (i) te be satisfied: 
2 
E, oe 38-58) Re 
ay ge es « 0.0246 (—3) (3518) 
(n, (&)) 10 


For the present set of experiments, (3.18) is satisfied for S1 es 2 oe 
whereas for the results of [1] it is satisfied for o1 Sj AUPE 
Since the present results are in agreement with those of [1], this 
limit may be extended to 51 =pOs225. Gon values ob ae 0.225 
the results are likely to involve an error larger than about 62. 
u,d 


The maximum value of cA in the present experiments 


was about 140, hence condition (ii) 4s sattisited: 
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Ce ae Wall Jet Region 

Considering the region where the flow is established to 
similarity, a first approximation is obtained if ,in spite of the wall 
shear, 1 the momentum is assumed to be preserved. That is to say 


the frictional force is small compared with the momentum flux involved. 
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With reference to Fig. 3.9, the momentum equation in the xy direction 


will give (neglecting ae 
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Thus, My is "preserved". On this basis, we could write: 
us a Y, OL P2Vox,) SE) 
But 
H i 
Se as 
xl ~ i Oe see leucs 
fe) fe) 
OL ert EO Geren Mea = const. Ps Hee Bue, by Eq.) 3 LOS 
oa 2b 
ne = Dime mp SOB on 
hence 
M = const.: QPU_-2b_ = const M 
al fo) 


Therefore: 


2 
u ,(0U,” 2b 5, ¥) 
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neglecting viscous effects. Hence: 


mi const (3.20) 
0 ¥x,/2b 
LEO 


Eq. 3.20 will hold approximately if the shear stress is 
relatively small. It is interesting to note that this equation is 


independent of H. Eq. 3.20 may also be written as: 
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ml Hw const 


(3221) 


re) re) Vx, /H 


ihe -data are plotted in this) form im Fie. 3.10, alone with the re— 


sults or [1]. Hq. .3.2land theretore Eq. 3.20 describe the data for 
x 


a certain distance - » beyond which the exponent of xX, appears to 
x 
differ from -+ . This value of - obviously depends upon the 


relative magnitude of ee From the above, we conclude that a more 


general statement of the problem is necessary, as follows: 
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3.3 Theoretical Development 


3.3.1 Free Jet Region 


We have seen in the previous chapter that (Eq. 2.14): 


“mL 2.40 
ue x 
ae 
10) 


To make this equation compatible with Eq. 3.4 which is a result of 


dimensional analysis, multiply with i nen: 


Pd es 2.40 
os as Vi 2 


H H/2b , 
; Drs ? ; aes , 
The quantity H/2b! is a virtual origin correction and generally 
fo) 
small in comparison with - . Thus, the above Equation may be re- 
written as: 
Yn fe 2.40 Cc, 
U 2b 
fe) fe) [xX 
H 


where Cy is an adjustment coefficient. The present data indicate 


a value of C) ~ 1.04, hence: 
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u 
a H a 2a (3.24) 


0 0 = 
H 


tioeiSevericied Incheon o.. and is, valid. for = <0. 65.55 che 


length scale for velocities is given by: 


The pressure Pa is very small in this region. The velocity field 
may be predicted by means of the similarity curve and the scales be 


Shave! el 
m 


Ooo 2 Impingement Region 
An aperomhara analysis can be carried out for the decay 
of u, as shown below. 
Assuming that uv is independent of x at small values of 


Wo (see t1e,0.2..20),.we could write’: 
duv 7 
(——) = k (3725) 


where “ke is a constant. But, 
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(using the continuity equation, a ae = = 0). Therefore: 
duv du 
ee Sten ae Ans 
eee Z (u 5. > 
y=o y=0 
since (vy) SMe 0 But, 
ENT saan 9 Seapine cnet cla 
ox on me eax re - 
hence, 
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u seer Dieta Git We 
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where ky is another constant. Using Eq. 3.4 we may write: 
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ik 2 Le x 
Gel oa 
oO (e) 
OL, 
oe 2 H DZ. x 
ne 2, = k, (ere are a2 (326) 


This’ is verified in Fie. 3.11; where the constants are found to be 


ke 25. 5s and ke = Oe. Hence: 
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(3.26A) 


Pigs equation tis plottedsin Fig. 2.2. “lt is seen that the unique 


function indicated in Eq. 3.4 may be approximated by two different 


x 


x 
laws. Between = 0.65 and aes O.7/5 2B SMOGOEM EreMISeOMm Eases, 


| 


However, neglect of this transition results in only small errors, 
For practical purposes, the boundary between the free jet and impinge- 
ment region may be defined as the intersection of Eqs. 3.24 and 3.26A, 


os 
O x 
which yselds 5-7 O./0%. 


H 
The maximum stagnation pressure on the wall, Po» may be 


predicted as follows. Applying the momentum theorem between the 


nozzle and the wall, we have (neglecting shear stresses): 
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* SAssimilareanalysis wae carried out for the case of circular amping— 
ing jets and was found to be in good agreement with the experimental 
results reported in [2] . 
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Then, 


p 
H 
ea 312 rt See Cm) 


pu,” /2 Oo K ow 
The constant oa is found (using Eq. 3.15) to be 0.143. The average 
value of Se by direct measurement was found to be 0.127. It appears 
that the integral tends to overestimate the value of ae hence it 
is suggested that oe be taken equal to 0.130. The maximum dis- 
crepancy would then be about 10%, This is attributed to the difiriculty 


of measuring wall pressure distributions precisely and in possible non- 


uniformity of the nozzle mouth. Using nase =ae0 19 Siig. a27. header. 
12 
—S- 5 = 7.7 (3.28) 
pu, /2 ro) 


To predict the axial pressure, Pm the variation of the 
total pressure, Poe on the axis ds- studied first... (“this is done, 
because if turbulence is reduced in the impingement region, then the 
total pressure will not vary excessively and with reference to Fig. 
S17, 2G mtent be possible towapproximate 1t by aystraienteliine.. noy— 
ing un? the variation of Pa could be found by subtraction. 


We have: 


To make this dimensionless, divide with Py : 
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P P pu 
T 
ee ee (3.29) 
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Dy Eq. 3.27. Substitutine in (3.29) we have: 


P Pp u 
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S Ps P fo) re) 


the data are plotted in this torm invtie. 3.125 where Pr 


is seen to vary mildly with x and the straight sline:s 


p 
Sr 0 (3.31) 
Pa | H 


appears to be a reasonable approximation. 


From Eq. 3.30 the axial pressure, Pa? is found as 


Ore by Las. 5.204, 3.27 "and 3925: 


ZZ 
a I (3°32) 


This is plotted in Fig. 3.3 and appears to be a reasonable description 
of the data for 7 »% 0.75. 

Equations 3.26 and) 3.32 cannot be eS eecsdere hold in a 
region very near the stagnation point, where the flow may approximate 
a potential pattern. However, the data indicate that such region is 
very small. 


With the information already available it is now possible 


to derive generalized equal-pressure contours. Let any contour value 


be n = 2. By Eqs. 2.12A and 3.11: 
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b Hb 2 
p P 
hence 
b ( 
y eet (2) y ES BS 5) 
H k H 
By GEG. 034.0% 
bp 
aio Bop (5) 9 
therefore 
ae) 
ate oe if 
H = Bop (5) ln 
In general: 
ae x 
A F(n, i (3535) 


Some contours’ are shown in Fig. 3.13. 
The variation of the maximum horizontal velocity component, 
La along the wall may be found by assuming the flow to be potential 


immediately outside the boundary layer. By the Bernoulli theorem: 
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Using Eqs. 3.28 and 3.15 and after some manipulations, we find: 


ey 72 
gqt8c5 (Ge) \ 


Teo e7 7 \ = (3.34) 


In Fig. 3.7, Equation 3.34 is seen to describe the data reasonably 
x 
well, for values of GZ & 0.4. Beyond this value, the decelerating 


effect of entrainment becomes appreciable and u decays steadily to 


ml 
assume wall jet variation further downstream. 
The wall shear stress, te: is found as follows. The 


equation of motion in the x-direction for plane flow is: 
Pe US ommreagichs -4 Bog 2a (3.35) 


where tT is the total shear stress constituting of both laminar and 
turbulent shear stresses, and wu, v, p are time average quantities. 


On the wall, this equation reduces to: 
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But oe = ae wee xy is referred to the coordinate 
x=H ik 


system x and denotes distance from the stagnation point along 


Tee ea) 


the wall. Hence: 
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eg = \ (ee (3.36) 
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The experimental data (Fig. 3.6) indicate that ee = 0. Using 


this observation we haves: 
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(3237) 
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The quantity Bo (1) is a constant, however, the available experi- 
mental data do not extend to the wall and it cannot be evaluated 


before hand. Let: 


' = 
85 (1) > 
Then, 
d C U - 2b 
( my Eg ae Peon oe 
dx K Z H H 
pw 
Substituting in (3.37) we find: 
C., ou,” Pe Dee n 
ee £(n) d 
ps z 5 a ( Hi ) f (n) dn 
pw 
O 
or, 
T : C a 
ee f sm dn 
2. 2b K 
pu, (210 pw 
chan 
nN 
uy. 9 H CG 
2(—) Sie Ae £(n) dn (3235) 
U 2b K 
Oo O pw 


oO 


But, 


HL - ty a _ 
“ Q 7 a Ssh ~ 


ith 1) : \ 
5 


78 


n n 2 ia] ie 2 
uf EQ) f aes! dj = ae fe ih d (kn) 
. £ vk A 
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oa 2 -Z 
ey a ee e dz . 
vk Vi 
O 
where z is a dummy variable. 
vn 9 
Lettines I.(vin) = oes e” dz we recognize I 
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to be the probability integral. Values of ro may be obtained from 


any tables of functions. Substituting in (3.38) we obtain: 
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we have finally: 
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The value of C, is evaluated from the shear stress data to 
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bewequal to O8l7. With “k = 0.693 and 


ogee 
q ne a 
we obtain: 
uy H = 
os 2, SO, yi y¥ I, (6.2 =) (3.40) 


Eq. 3.40 is) plotted in Fig. 3.8 where it is seen to describe the data 
% 

reasonably well, for values of a Jess than 074. Beyond this poine 

the normal pressure gradient cannot be predicted by similarity con- 


siderations. From (3.40) we find: 
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eee oC) AD 


The value 0.17 for Ci indicates that oS = 0.053. That is) to-say 
the axial pressure gradient at the wall is very small, owing to 
occurrence of nearly potential flow at the neighborhood of the stag- 


nation point. 


Sooo eNall Jet Reelon 


After a certain distance downstream of impingement,the flow 
becomes established, so that the velocity distributions at different 


x, stations may be reduced to a single similarity curve (Fig. 2.30): 
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where by is the value of yy at which: 


du 
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u aE ual and cs <o Om, 
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The Reynolds equations, after the usual boundary layer approximations 


read (plane turbulent flow): 
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du ov 
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i 1 
ee 8 0 ’ ' ' 4 , 
where es pu; Vy and uj » Vy are the fluctuating components 
of velocity. 
(ie We an D are time average quantities. 
i a ih 
The pressure Py is very nearly atmospheric throughout the 
Py 
Jet, hence eae Q . Therefore: 
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Due to similarity; i = ud £(n,) . Assuming that the shear stress 


profiles may be made similar by means of the maximum dynamic pressure, 
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we could write: 
= 2 


Substituting the similarity expressions in Eqs. 3.43 and 3.44 we 


obtain after some manipulations: 
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where Gi> G are functions of n 


9 only. Since the variables x 
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Ny are formally independent: 
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Eq. 3.46 yields: 
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where C is a constant and x is measured from a virtual origin 
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for the length scale by : 


The linearity of by is verified in Fig. 2.33. The value 
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of C. 160.04. ” Thus: 


Ch = CC aD? R ) Jand C, = Cy ( mo RY 
The present results suggest that Ch does not depend on ie peels 
) 
is also supported by the results of [1] . However, the value of C. 
from (| llyis 0.10. which andicates: that 
dc. 
= — > ‘ 
C. C, (RA) and aR. 0 (3.49) 
Using Eq. 3.47, we find: 
= k x,* (3.50) 


where k, a are constants. 


Returning to Eq. 3.44 we rewrite (using Eq. 3.43): 
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rule, we have: 
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Due to similarity the value of the integral is a constant, say F. 


Hence: 


Usine (3.48) and (3.50) we finds 
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where Ce is a local skin friction coefficient, defined as: 


” _ 
- 7 | re | , a oy - a 
Tam ‘eral i235! wit) Has rc = 
‘ 0 ™ f A ae4. Wola neil A 


at 
a: 


a 


84 


Eq. 3.51 indicates that the local skin friction coefficient 
u x 


: ° * 
is independent of Xi Figure 3.14 shows aaa plotted against ee 
ml 
xa u 
After a certain value of ne ete becomes indeed a constant. This 
ml 


is in agreement with the conclusions of Schwartz and Cosart [18]. 


Sigalla (19] reports a dependence of Ce on the local Reynolds number 


way 
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however, this variation is very "mild" so that in substance 


Chere as no contradiction. § The constancy of | C in the wall jet may 
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also be derived in a different way, which can indicate its dependence 
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on other parameters such as RG and oD If the defect law is 
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assumed to be valid in the boundary layer for some range of Yy> then: 
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Or. 
ee ree ae ) 
— = _—_— ia] 
“ml “ml : : 
on, 
u, £(n,) => ih 
Us 
This suggests that ee and therefore Ces is independent of Xy> 
ml 
and depends only on the form of the functions f and hy . Functions 


£(n,) and h,(,) are unlikely te change with the dimensionless 


height, however, it is possible that they change slightly with 


i 
2D 
fe) 
Reynolds number. 


Values of Ce are plotted against R, ag aes US Aree ali 2 

is seen that Ce decreases with Ro whereas any dependence on iia 
fo) 

appears to be very small. To evaluate fully the variation of Ce 


with Ro more experiments may be required. As a conclusion we may 
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state that un varies as ua with the constant of proportionality 


being deperident on R, x 
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suggests that as long as Ce Pe AW ees -4 . We have seen already 
dc, dc. 
that = = — —— 
a C. CR)» Ce C.(R,) and aR SOM, aR, <hcOie We 


may therefore conelude that: 
a d|a| 
|a| = a(R.) and aR < 0 (3303) 
Oo 
Eq. 3.52 should be used only as an approximation for values of the 
exponent a. The actual value of F should be generally smaller 
than 0.70 due to negative velocities occuring outside the jet. 


The following table shows values of the exponent a com- 


puted from Eq... 4.52, using FF = 0270. 


TABLE IT EXPONENT a 


Source 


Present 


Experiments 
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The present results were evaluated using a = 0.60 in 
Figures 3.16 and 3.17. The results of [1] are plotted in the same 
manner in Fig. 3.18 using a = 0.52. The exponential variation is 


verified for values of xy greater than about 70 nozzle widths. 


From the above, an equation for ul is derived: 
u x 
milo ik a 
i C Gp Co) (ae 5). 
Oo fe) 


Recalling Eq. 3.22, we conclude that in general: 


We have seen already that a = a(Ro)- Values of C are plotted 


against oles in Fig. 3.19 where a single line may be drawn through 


2b 5 
the data. Although these values of C are tentative it appears that 
C depends little on Ro » the coetticient Cy fs ~a virtual origin 
H 
correction and appears to depend on both Ro and > For the 
fo) 


present experiments C_* 35, whereas for the results of [1] Cy = 0. 
A weakness of this analysis lies in the fact that the 
virtual origins for velocity and length scale do not coincide. This 
is very common in jet flows and at present there is no plausible 
explanation for the discrepancy. 
In spite of this the analysis may be retained to serve as a 


guide for evaluation of experimental results. Moreover, at Reynolds 
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numbers higher than about 4 x 10° the distance between the virtual 
origins becomes so small that it could be neglected altogether with- 
out causing appreciable errors in computing oe and by 5 


With the results already available it is possible to draw 


a general outline of the impinging jet. This is done in Fig. 3.13. 
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CHAPTER IV 
SUMMARY AND CONCLUSIONS 


An experimental and analytical investigation was carried 
out regarding the mean flow characteristics of plane turbulent jets 
impinging Me smooth walls. 

In agreement with previous investigations three distinct 
flow regions were found to exist. These are: 

(i) The Free Jet region. 
(ii) The Impingement region. 
(iii) The Wall Jet region. 
The properties of region (i) were found to be almost identical with 
the corresponding properties of the free jet. Thus, upstream of a 


DOD Ca Xo it was found: 


E 
N 
is 
=) 


b 
Me Re Ue x e e 
= ete O.7H and a 0.100 + 0.15),x, 0.65H 


The static pressure in excess of the hydrostatic was not zero, but it 

was found to be generally small. The velocity profiles were found to 
u : 

be similar in the sense that ae = Sr) . The function f£ may be 


m u 
approximated by an error distribution. It is also in good agreement 
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with Tollmien's classical solution. 

Flow properties in region (ii) were investigated by means 
of a Pitot cylinder which was calibrated first. The calibration was 
based on the assumption that within a range of tube Reynolds numbers 
the pressure distribution around the cylinder does not change appreci- 
ably with Reynolds number. This was verified within the range 
Ro = 1400 to 5500 for values of polar angle $ from O to 70°. 
Caution should be excercised if it is desired to use the calibration 
charts presented here at Ro values far below or beyond the above 
range. In view of the fact that the pressure distribution varies 
appreciably with Ro in the wake zone, it is suggested that for 
RL >> 5500 (or Re <<) £4005) the present calibration charts be used 


Hg 


with the restriction || ¢ 30°. A lower limit of R, Should be in 
the order of 300. 

In the impingement region the vertical velocity distributions 
were TOund sto remain Similar for x74 .0.65Hy For (x 2.0.65 the spro— 
files are not similar although deviation from the similarity curve is 
not extreme. For crude approximations, the similarity curve could 


be used. The length scale, be was studied experimentally. The 


pressure profiles were found to be similar in this region, so that: 


The length scale, Dae was studied experimentally. 
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The velocity and pressure scales are given by: 


u 

_m aH We x 

7 Rc Sg ee seh 
Oo Oo 

p 

ae oe Slee tel x > 0.73H 
Ss 


The total pressure, Pps on the axis is given by: 


The variation of u, was obtained on the basis of the 
experimental observation that the flux of x-momentum in the y 
direction is almost independent of x near the axis. This physical 
concept was applied to the circular jet yielding good agreement with 
available experimental data. The above equations are likely to be 
erroneous in a small region very near the stagnation point, because 
the flow tends to become potential in this region. 


The stagnation pressure, Ps? may be predicted from: 


Ps H 
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pu, 72 O 


The wall pressure distribution was found to agree well with the 
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function: 


The maximum horizontal velocity is given by: 


u 2 
1 A = 

ao = 2.77 f 1 -} for! 4. 510.44 

U 2b Pp 1 

) fe) s 
The wall shear stress is given by: 

T x 

ee ee ete) eag Fe oO 
2 2b P H al ; 

pu, f2 Oo 


x x 
where bo? =) is the probability integral between O and ee : 
According to Patel's criteria for favourable pressure gradient flows, 


the experimental results for wall shear stresses are likely to in- 


volve errors larger than 6% for xy aS WAS )ele e 
In region (iii) the flow is established to self-similarity 
pal ah 
in the sense: — = f(n,), nn. ==. after approximately 70 
wal ih ul by 


nozzle widths from the stagnation point. 


The velocity scale is given by: 
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where a is a function of Reynolds number, Ro » and is generally 


a 
smaller than - 7° As Ro increases, a tends to ee 


7 and 16 


could be taken as - $ f£oV Ro Pan) 000. 


C is a parameter depending generally on a and R, 
fe) 


however, the dependence on Ro was indicated to be weak. 
Ce is a virtual origin correction. Generally, S 


depends on ae and Ro . For the present set of results = 35 
fo) 
4 


whereas for the results of [1], where Ro So 4 XielO) £s. ~ 0. 


The length scale, by o8, Us, piven bys 


by = COs a) 


where x is a virtual origin correction and Ch was found to depend 


only on Ry : 


The shear stress was found to be proportional to the square 


of the velocity scale, ua? ie: 
ig x 
ee Oe oe Oss 
219 f£ H 
Punt 


These local skin £fietion coeriicient, Ce ailSsa 1unceion, oF Ro only 
dc 


ey 5 
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and 
The results reported in this thesis should be applicable for 
impingement heights not smaller than the length of Mtheipotential core 


Of the Snutial jet. 
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